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The recent Boomerang and MAXIMA data on the cosmic microwave background suggest a large 
value for the baryonic matter density of the universe, ujh ~ 0.03. This density is larger than 
allowed by standard big bang nucleosynthesis theory and observations on the abundances of the light 
elements. We explore here the possibility of accommodating this high density in inhomogeneous big 
bang nucleosynthesis (IBBN). It turns out that in IBBN the observed D and Yp values are quite 
consistent with this high density. However, IBBN is not able to reduce the ^Li yield by more than 
about a factor of two. For IBBN to be the solution, one has to accept that the ^Li plateau in 
population II halo stars is depleted from the primordial abundance by at least a factor of two. 
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I. INTRODUCTION 

The most accurate way to estimate the average den- 
sity of baryonic matter in the universe, ph, has for a 
long time been big bang nucleosynthesis. In standard 
big bang nucleosynthesis (SBBN) ^ the calculated pri- 
mordial yields of the light elements depend only on the 
baryon-to-photon ratio r] = Ub /rij. Future precision mea- 
surements of the fluctuations in the cosmic microwave 
background (CMB) will provide another way of measur- 
ing the baryon density. The shape of the angular power 
spectrum of these fluctuations will depend on a number 
of cosmological parameters, among which is the baryon- 
to-photon ratio. In this context it is usually given as the 
baryonic contribution to the critical density times the 
Hubble constant squared, uib = ^th^- 

These two measurements rely on completely different 
physics, and if they agree it is an important confirmation 
that we have the right understanding of the early uni- 
verse. Assuming that there was no significant entropy 
production in the universe between BBN and recom- 
bination, the two parameters cob and rj are related via the 
present temperature of the CMB, T = 2.725 K, by 
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Two recent balloon-borne experiments. Boomerang 
1^ and Maxima- 1 [Q, have now provided us with the 
first measurements of the CMB angular power spectrum 
which are of sufRcient quality that an estimate of ujb can 
be made from them. The result cjh ^ 0.030, is sig- 
nificantly higher than the SBBN result. Indeed, SBBN 
clearly cannot accommodate as high a baryon density as 
LUb ^ 0.030 §. 

The SBBN yields for cut = 0.030 are Yp = 0.251, 
D/H = 1.7 X 10-5, 3He/H = 8.4 x 10-^ and ^Li/H = 
7.6 X IQ-^". With the exception of ^Li, the uncertainty 
in these SBBN yields is much less than the uncertainty 
in the primordial abundances derived from observations. 



There is much debate about chemical evolution and 
systematic effects in observations. For the primordial 
helium abundance there are two competing estimates, 
the "low "He" 0, = 0.234 ± 0.003 and the "high 
''He" |], = 0.244 ± 0.002. The difference is largely 
due to the method of estimating the present abundances 
from the observed line intensities, suggesting that the 
systematic uncertainty may be larger than the ±0.005 
usually assumed 

Buries and Tytler claim to have established the 
primordial deuterium abundance as D/H = (3.3±0.25) x 
10~^ ("low D") based on Lyman-series absorption by 
high-redshift clouds. This is based on a detected low deu- 
terium abundance in three such systems and upper lim- 
its from others. There remains one such system, where a 
high deuterium abundance, D/H 2 x 10"'^ ("high D"), 
has been observed It may be that the accuracy of 
such observations has been overestimated We re- 

fer the reader to recent reviews |ll3| for further discus- 
sion and adopt the observational constraints Yp = 0.228- 
0.248, D/H = 2.9-4.0 x lO^^ These lead to the 

SBBN range cJb = 0.004-0.021 from Yp and cob = 0.018- 
0.022 from D/H. The "low D" of Buries and Tytler gives 
ivb = 0.020 ±0.01. 

The "Spite plateau" ||l^ of ^Li abundance in popula- 
tion II halo stars provides us with the best estimate of 
the primordial lithium abundance. Bonifacio and Molaro 
[ p^ obtain a present lithium abundance log]^o(^Li/H) = 
—9.80 ± 0.012 ± 0.05 for these stars, and argue against 
any significant depletion from the primordial abundance, 
based on the lack of dispersion in the data. Pinson- 
neault et al. |17| estimate a depletion factor of 0.2-0.4 
dex from rotationally induced mixing. In a more re- 
cent study Ryan et al. [ p^ obtain a mean abundance 
logiQ(''Li/H) — —9.88 for the Spite plateau, and argue 
that the narrow spread, less than 0.02 dex, limits deple- 
tion by rotationally induced mixing to less than 0.1 dex. 
Moreover, they observe a slight trend with metallicity, 
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suggesting a galactic contribution, so that the primor- 
dial abundance could be lower than the observed abun- 
dance. Based on this, Ryan et al. derive the range 
logioC^Li/H) ^ -10.04. . . -9.72 for the primordial abun- 
dance. This very tight upper limit to primordial lithium 
is in conflict with the "low D" estimate jl^], ujb = 0.020, 
which corresponds to log]^Q(^Li/H) = —9.5 in SBBN. 

At present there seems to be a consensus p7| , pO| , ^ , 
that the primordial abundance cannot be larger than 
the observed value by more than about 0.4 dex, so that 
models producing a primordial lithium abundance above 
logio(^Li/H) = -9.4, or ^Li/H = 4 x 10'^°, would be 
ruled out. This gives an upper limit cub < 0.022 in SBBN. 

While the upper limit to ujb from Yp is rather soft, the 
deuterium and lithium abundances are in clear conflict 
with uJb = 0.03. 

This CMB measurement of ujb is rather preliminary 
and not very accurate, and thus it would be premature 
to discard the SBBN estimate. Since the CMB measure- 
ment is a simultaneous fit of many cosmological param- 
eters, the lack of prominence of the second peak in the 
power spectrum could have another explanation. Sug- 
gestions include a significant tilt p^ , an isocurvature 
component p^ , or a feature in the primordial per- 
turbation spectrum, and topological defects p^ ]. 

Future satellite experiments, MAP |2^ and Planck 
, will give us a better CMB determination of Ub- If 
the situation persists and the CMB value of uJb turns out 
to be LOb ^ 0.03 or higher, one has to consider abandon- 
ing SBBN in favor of non-standard BBN There are 
many modifications suggested to standard BBN, some of 
which can accommodate a higher baryon density |p9| . 

We consider here one such possible modification, in- 
homogeneous big bang nucleosynthesis (IBBN) . In IBBN 
one assumes that the baryon density was inhomogeneous 
during BBN. This inhomogeneity could be an initial con- 
dition resulting from unknown physics in the early uni- 
verse, or it could be the result of some physical event, 
like a phase transition, occurring before BBN. 

The nature of IBBN depends on the distance scale 
of the inhomogeneity. Especially interesting is the case, 
where this distance scale is of the order of the neutron 
diffusion length during BBN |^^. Then neutron diffu- 
sion will lead to an inhomogeneity in the neutron-to- 
proton ratio. If the density contrast between high- and 
low-density regions is large enough, this ratio can exceed 
unity in the low-density region. In this case IBBN leads 
to a decrease in the ^He yield and an increase in the D 
yield, favoring larger uJb- However, the ^Li yield is more 
likely to increase than decrease from the SBBN result for 
the same average uJb- 

In this paper we focus on the question, whether the 
high value LUb = 0.030, corresponding to ?7io — 8.22, 
would be acceptable in IBBN. There have been many pa- 
pers on IBBN [Q. However, because of the large number 
of parameters in the IBBN model, the parameter space 
has not been covered thoroughly in them. Focusing on 
one particular value of ijio we reduce the dimension of 



the parameter space, allowing us to search the space of 
the remaining parameters more completely. We can thus 
give a more specific answer to the above question. 

II. RESULTS 

We have updated our IBBN code |^l[ with the latest 
compilation of reaction rates |^2| . The new reaction rates 
have little impact on the isotope yields. The largest effect 
is on ^Li production, which is reduced by 5%. Vangioni- 
Flam et al. estimate the uncertainty in ^Li yield due 
to reaction rates to be 30%. The actual uncertainty may 
be even larger, since their method of error estimation 
allows error compensation. 

We assume a simple initial density profile, so that there 
is a high density region, with one baryon density (r?high), 
and a low-density region with another baryon density 
(jliow)- Diffusion will then smooth out this initial den- 
sity profile. We assume a spherically symmetric geome- 
try, allowing us to study two different geometric config- 
urations: 1) centrally condensed (CC), where we have a 
spherical high density region surrounded by a low den- 
sity region, and 2) spherical shell (SS), where we have 
a spherical shell of high baryon density inside which we 
have low baryon density. The first case approximates a 
situation where the high density regions have a rather 
compact shape, with a low surface-to-volume ratio, the 
second case a situation where the high density regions 
have a more planar geometry, with a large surface-to- 
volume ratio. 

Once the average baryon density rj is fixed, we have 
three remaining parameters, the density contrast R = 
77high/?7iow, the volume fraction fy of the high-density re- 
gion, and the distance scale r. There is an "optimum 
value" Topt of the distance scale, related to the neutron 
diffusion scale, where the effect of the neutron excess is 
maximized. For very high density contrasts, fyR ^ 100 
the results become independent of R. For a more detailed 
discussion, see Ref. Q. 

It has been customary in IBBN work to choose R suffi- 
ciently large to get rid of the i?-dependence. This reduces 
the dimension of the parameter space by one, and maxi- 
mizes the effect of the inhomogeneity. We first follow this 
practice by fixing fyR = 100, and show the IBBN yields 
of "^He, D, and ''Li in the remaining two dimensional pa- 
rameter space (/«,r) in Fig. |^ for the CC geometry, and 
in Fig. ^ for the SS geometry. 

All our distances given in km refer to comoving dis- 
tance at T = 1 MeV. Note that 1 km at 1 MeV corre- 
sponds to 8.5 m at 100 MeV and to 1.9 x 10^"* pc today. 

We see that IBBN leads to acceptable ^Hc and D yields 
for UJb = 0.030. This occurs for distance scales which are 
near the "optimum distance scale" for IBBN. Moreover, 
choosing the "high" observed value for Yp, and the "low" 
observed value for D/H gives us agreement for moderate 
(not very small) values of /„, where the ''Li yield is the 
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lowest. Choosing a lower Yp and a higher D/H would 
move the agreement towards smaller /„, where the ^Li 
yield would be significantly larger than in SBBN. How- 
ever, for the SS geometry, somewhat higher D/H values 
would be fine. 



For SBBN with ujb = 0.030, all three observables, Y, 



D/H, and ^Li/H are outside their observed ranges. While 
the inhomogeneity at the right distance scale had the 
effect of moving Yp and D/H into the observed range, 
this is not true for ^Li/H. For the CC geometry IBBN 
does not reduce the ^Li yield appreciably. For the SS 
geometry there is a reduction in ''Li/H for not too small 
fv, which could arguably be sufficient. 

To see what the effect of the remaining parameter, 
R has on this lithium problem, we choose to track the 
"lithium valley" of Figures and |. In ^ we derived 
the dependence of the optimum scale Topt on /„ and rj, 



,1/3 

-2/3 



r-opt 



(CC) 
(SS). 



(2) 



From Figures and I we see that the minimum 'Li yield 
occurs for slightly smaller distance scales than the mini- 
mum ^He yield. We choose the proportionality constant 
so that r(/^) tracks the low ^Li yields. Thus we have 
taken 



r = 25km 



1/3 



(1 - fv) 



(CC) 



(3) 



= 32.5km- 
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(SS), 



(the dotted lines in Figures 1 and 2). Wc now allow R to 
vary independently of fy and show the results in Figures 
I and I 

To check that the lithium valley does not shift in r as 
the density contrast R is reduced, we repeated this for 
nearby distance scales. 

These figures quantify the effect of the density contrast 
R in IBBN. We see how the effect indeed saturates for 
fyR > 100. For IBBN to have the desired effect of bring- 
ing Yp and D/H into agreement with observations for a 
high baryon density (here Wf, = 0.03), we need a density 
contrast of R ^ 4//,j. Varying R does not have much 
effect on the lithium problem. The smallest ''Li yield we 
get with acceptable Yp and D/H is ^Li/H = 7.4 x 10" i° 
for the CC geometry and ^Li/H = 4.0 x IQ-^" for the SS 
geometry. 

We found that IBBN with the CC geometry is not able 
to reduce the ""Li yield significantly, but the SS geometry 
can reduce it by about a factor of two. Thus one could 
argue for a marginal IBBN case, where we have an ini- 
tial inhomogeneity approximating the SS geometry, i.e., 
the surface-to-volume ratio of the high-density regions 
is relatively large. With a high-density volume fraction 



fv = 0.2-0.5, density contrast R > 30, and a distance 
scale r — 35-70 km at 1 MeV, we get Yp near the upper 
limit 0.248, D/H comfortably in the observed primordial 
range, and ''Li/H ^ 4 x 10""'^°. With the uncertainties 
in the ^Li/H observations and BBN reaction rates lead- 
ing to ^Li this is in agreement with observations if one is 
willing to accept a depletion factor of about two. 

We are, however, not aware of a well-motivated sce- 
nario for producing a baryon inhomogeneity both with 
this geometry and this distance scale in the early uni- 
verse. 

A baryon inhomogeneity could be generated during in- 
flation as an isocurvature perturbation, or after inflation 
in some phase transition. In the latter case, the distance 
scale must be less than the horizon size during the phase 
transition. At the electroweak phase transition, the hori- 
zon size is about 3 km, and the expected distance scale is 
significantly smaller than this, too small to have a large 
effect on BBN |l|- 

The QCD transition is the most studied candidate for 
generating the inhomogeneity for IBBN |Q. The hori- 
zon size at the QCD transition is about 2000 km. In a 
first-order transition, the resulting baryon density con- 
trast could well be large enough, and the distribution of 
distance scales would be fairly narrow. In homo- 
geneous thermal nucleation, the distance scale depends 
on the latent heat L and the surface tension a of the 
phase transition. Typical estimates from lattice QCD 
calculations give distance scales which are much too 
small for our scenario, although combinations of L and 
a leading to sufficiently large distance scales may not be 
completely ruled out [Q. Heterogeneous nucleation due 
to impurities could also lead to larger distance scales, of 
about 1 km |^^. (Note that in the literature these dis- 
tances are often quoted as comoving at 100 MeV, and 
thus smaller by a factor of 1/117). 

Because of primordial perturbations the nucleation of 
the phase transition does not take place in homogeneous 
conditions. The cosmic fluid is undergoing acoustic os- 
cillations at subhorizon scales. Since the nucleation rate 
is extremely sensitive to the temperature, thermal nucle- 
ation would occur only at cold spots. This could lead to 
a distance scale of the order of 1 km , which is large 
enough to affect nucleosynthesis. However, the inhomo- 
geneities produced by the QCD transition appear likely 
to resemble the CC geometry more than the SS geometry. 



III. CONCLUSIONS 

It is intriguing, that the "high" baryon density Wf, ~ 
0.03, suggested by recent CMB results fits so well with 
the "low" D/H and "high" Yp observations in the IBBN 
scenario. However, the ^Li yield is significantly higher 
than the standard estimates for its primordial value. 

For the SS geometry, IBBN can reduce the lithium 
yield by about a factor of two in the best case, but for 
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ujb = 0.03, we still get ''Li/H = 4x 10~^°, requiring deple- 
tion by at least a factor of two in the Spite plateau. This 
may be marginally acceptable, especially if one allows for 
the uncertainty in the ^Li yield due to the reaction rates. 

We note for comparison, that the widely accepted 
"low-D" determination ujb = 0.020 in SBBN yields 
"^Li/H = 3.4 X 10"^° and thus also requires significant 
lithium depletion p9| . 

The better developed scenarios for producing the kind 
of inhomogeneity needed for IBBN tend to lead to the CC 
type of geometry. For this geometry, IBBN is not able 
to reduce the lithium abundance significantly below the 
SBBN result, so that depletion by a factor of four would 
be needed. This seems to be considered unacceptable at 
present. 
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FIG. 1. The light element yields (a) Yp, (b) logi(,(D/H), 
and (c) log]^Q('^Li/H) as a function of the high-density volume 
fraction fv and the distance scale r, for the CC (centrally 
condensed) geometry and ut = 0.030. The density contrast 
R is fixed hy f^R= 100. The observational Umits Yp < 0.248 
and D/H = 2.9-4.0 x 10~^ are indicated by thick lines (both 
are shown superposed on the ^Li yield in (c)). The dotted 
line shows the r(/„) relation, "optimum valley", chosen for 
Fig. 1^. The limit of small r or large /i, gives the SBBN result. 



5 



0.01 , 0.1 1 



V 




0.01 , 0.1 1 



FIG. 2. Same as Fig. |l|, but for the SS (spherical shell) 
geometry. The dotted line shows the r(/„) relation, "optimum 
valley", chosen for Fig. 0. 
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FIG. 3. The lithium yield, logj^Q(^Li/H), as a function of 
the high-density volume fraction /„ and the density contrast 
R for the CC geometry. The distance scale r corresponds to 
the dotted line in Fig. hi. The observational limits Yp < 0.248 
(above the solid line) and D/H = 2.9-4.0 x lO"'^ (between 
the dashed lines) are indicated by thick lines. The dotted 
line corresponds to the relation fvR = 100 used for Fig. ^. 
Taking i? = 1 or /„ = 1 gives the SBBN result. To bring out 
the slight decrease in the^Li/H yield, we have included the 
logjg(^Li/H) — —9.13 contour in this figure. 
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FIG. 4. Same as Fig. |[ but for the SS geometry. 
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